Edge magnetism of Heisenberg model on honeycomb lattice 
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In our previous study, a single-branch of ferromagnetic magnon with linear dispersion is shown 
to exist near the (uncompensated) zigzag edge for Heisenberg model on honeycomb lattice. Here 
we develop a field-theory description for the edge magnon and find its dynamics is captured by the 
one-dimensional relativistic Klein-Gordon equation. It is intriguing that the boundary field theory 
for the edge magnon is tied up with its bulk counterpart, described by the two-dimensional Klein- 
Gordon equation. Furthermore, we also reveal how the parity symmetry relates evanescent modes on 
opposite edges in a honeycomb nanoribbon. By employing alternative methods, including Schwinger 
bosons and density-matrix renormalization group, we also demonstrate that the relativistic edge 
magnon is robust even when the Neel order in the bulk is destroyed by quantum fluctuations. 
The edge magnon is a direct consequence of uncompensated edge and may be verified in realistic 
materials by experimental probes. 



Introduction. Exchange interactions between lo- 
cal magnetic moments, often described by Heisenberg 
model[l] and its derivatives, lead to rich and some- 
times exotic phases in quantum magnetism[2]. For in- 
stance, the excitation gap in the integer-spin chain pro- 
posed by Haldane[5 stimulates theoretical investigations 
and is later verified in experiments [H [5]. Antiferromag- 
netism in two-dimensional square lattice has been stud- 
ied extensively because of its adjacency to unconven- 
tional superconductivity in cuprates[5] and iron-based 
materials [5. Recent breakthrough shows that 5 = 1/2 
Heisenberg model on Kagome lattice exhibits exotic spin- 
liquid ground state[8] due to strong quantum frustra- 
tions. Moreover, it has been demonstrated that superex- 
change interactions between ultracold atoms can be re- 
alized in optical lattices |^. It may provide a different 
route to understand various ground states of the Heisen- 
berg model on different lattice structures. 

It is known that boundary effects give rise to frac- 
tionalized excitations in integer-spin chains |10l lllj but 
are less studied for spin systems in higher dimensions. 
The importance to understand the boundary effects in 
Heisenberg model is echoed by potential edge mag- 
netism in graphene nanoribbon |12f[l8] . Monte Carlo 
simulations [TS] demonstrate mean-field like ferromag- 
netic moment near the zigzag edge of graphene nanorib- 
bon. Unlike the usual ferromagnetic magnons with 
quadratic dispersion, these sharp spin excitations near 
zigzag edge shows linear dispersion. Recent spin-wave 
calculations [20^ for the Heisenberg model on honeycomb 
lattice show that the dispersion of the ferromagnetic edge 
magnon is indeed linear. 

The emergence of edge magnetism can be understood 
by Lieb's ferrimagnetism. However, it remains puzzling 
how the ferromagnetic edge magnon acquires linear dis- 
persion, which is clearly related to the antiferromagnetic 
correlations in the bulk. In fact, the absence of quadratic 



FIG. 1: Neel state in a honeycomb nanoribbon with zigzag 
edges. Spin orientations (purple arrows) on sublattice A 
(green dots) and sublattice B (red dots) are opposite to each 
other. The choice of unit cell is highlighted by the shaded 
blue circle. 



dispersion hints that the boundary field theory for the 
edge magnon must tie up with the bulk, just like the 
chiral states in quantum Hall effects. In this Letter, 
we study edge magnetism of the Heisenberg model on 
honeycomb nanoribbon with zigzag edges as shown in 
Fig. 1. Within spin- wave approximation, magnons in 
the bulk follow the two-dimensional Klein-Gordon (K- 
G) equation, characterized by the spin-wave velocity Cf, 
and the effective mass mj,. The boundary condition 
gives rise to evanescent modes and imposes constraint 
on the propagation along the edge and the (imaginary) 
momentum in the transverse direction. It is quite re- 
markable that the edge magnon is described by the one- 
dimensional K-G equation with parameters Cg and m,, 
related to their bulk values. Our derivations provide nat- 
ural explanation for the linear dispersion and reveal the 
connection between the boundary field theory for edge 
magnons and its bulk counterpart. Though the above 
conclusions are drawn from spin-wave theory, we further 
extend the calculations by Schwinger boson approach. 



where the bulk no longer has long-ranged Neel order, 
and find the edge magnon exhibit qualitatively the same 
behavior. Meanwhile, we also perform density-matrix 
renormalization group (DMRG) calculations to demon- 
strate that the spin-singlet ground state in nanoribbon is 
actually very close to the Neel state. 

Harper equations. To explore the boundary effects for 
Heisenberg model on honeycomb nanoribbon, we first 
write down the Hamiltonian for exchange interactions, 

[r.r') a 

where (r, r') is denoted the nearest-neighbor pairs and 
Jx — Jy ~ J, Jz = iJ are exchange couplings with 
anisotropy 7 > 1. The spin operators can be repre- 
sented by Holstein-Primakov (HP) bosons. Assuming the 
ground state has Neel order as shown in Fig. [T] the in- 
teractions between these bosons can be ignored and the 
effective Hamiltonian within spin-wave approximation is 

Hsw = JSJ2 f^ {bAir)bAir) + bUr')bB{r')) 



{r.r') 



+b\{r)bl{r') + bB{r')bA{r) 



(1) 



where S is the magnitude of the spin and bj^/g are an- 
nihilation operators for HP bosons on sublattices A/B. 
Since the spin- wave Hamiltonian is bilinear in HP bosons, 
it is more convenient to write down its equivalent equa- 
tions of motion in first-quantization language. Following 
the same steps developed in Ref. 20J, the dynamics is 
described by the coupled Harper equations, 

XipA{r) + QY,^*B(r + S,)^idtipA{r), (2) 

Si 

QY,^A{r + S,) + Xip*Bir) = -idt^pUr), (3) 



where r denotes the lattice sites for honeycomb lattice 
and Si are the vectors pointing to the nearest neighbors. 
The wave functions on different sublattices are (/3a,Vsi 
where the conjugation arises from opposite spin orienta- 
tion. The parameters are Q = JS and A — zjQ, where 
z is the number of nearest neighbors. It is important 
to emphasize that the magnon carries quantum number 
ASz = =f1 and sets the normalization condition, 



AS.^-Y,[\ipA{r)f-\ipBir)\')^Tl 



(4) 



The above Harper equations can be solved exactly, de- 
livering a single-branch ferromagnetic magnon near the 
zigzag edge with linear dispersion. In the following, we 
would like to develop general field-theory descriptions to 
explicitly reveal the connection between magnons in the 
bulk and those at the edge. 



Field theory for the bulk. In the field-theory limit, 
we introduce the smooth-varying fields, 4>{^(r,t) = 
V ^\k\<Aa "^aC^'O e*'"'', where the momentum summa- 
tion is restricted to the vicinity of A; = with a cut- 
off Ac. For these smooth- varying fields, spatial variable 
r = (x, y) can be treated as continuous and no longer 
restricted to the lattice sites. As a consequence, spatial 
derivatives are well-defined. Making use of the displace- 
ment operator, e°''^4>{r) = (j){r + a), the Harper equa- 
tions can be represented in the matrix form, 



A 
Qh{dx,-dy) 



Qh{d^,dy 
A 



idt 



(5) 



where h{d^,dy) = 2e''^''^y/^ cosh{ad^/2) + e'^^'^y^^ for 
honeycomb lattice. Keeping the lowest order in gradient 
expansions and eliminating the field (j)*g, the dynamical 
equation solely for the field (f>A can be derived. It is 
not surprising that the effective field theory turns out to 
be the well-known Klein-Gordon equation in two dimen- 
sions. 









dy"^ 



2 2 



0^(x,y,i) = O. (6) 



The spin-wave velocity in the bulk is c;, = ■\/3/2 Qa and 
the effective mass is m^ = 2-^/7^ — l/(Qa^). One can also 
eliminate the field (pA and show that (jfg also satisfies the 
same K-G equation. These results are not surprising be- 
cause antiferromagnct in low-energy limit is relativistic. 
Without the annoying spin kinematics, spin operators 
can be viewed as canonical bosons and K-G equation is 
a natural description for relativistic bosons. 

It is important to keep in mind that (J)a and (jf^ are 
indeed antiparticles to each other. As required by rela- 
tivity, they always appear in pairs and explain the double 
degeneracy for magnons in an antiferromagnct. Further- 
more, when the anisotropy disappears, 7 = 1, the excita- 
tion gap for the magnon mi,c\ also disappears as expected 
from Goldstone's theorem. 

Field theory for the edge. One can also introduce the 
smooth-varying fields on the edge and applies the same 
techniques to derive the boundary field theory. Since 
our goal is to demonstrate the connection between the 
field theories in the bulk and on the edge, it is wise to 
write down the field-theory presentation for the boundary 
conditions. For the honeycomb nanoribbon considered 
here, at the upper edge {y = Ly) where the outmost sites 
belong to sublattice A, the boundary condition gives the 
constraint 7(^^(0;, Ly)-|-(/?^(a;,Ly-|-aj^) =0. On the other 
hand, for the lower edge {y — ~Ly), the outmost sites 
belong to sublattice B and the boundary condition leads 



to -y^six, -Ly) + (p*a{x, -Ly 



0. As long as the 



transverse width Ly is finite, edge magnons on opposite 
edges entangle together and complicate the problem. For 
simplicity, let us temporarily assume that the transverse 




FIG. 2: Dispersions for edge and bulk magnons, (a) in 
isotropic limit ^ — 1 and (b) with slight anisotropy 7 = 1.01. 
The dispersion of edge magnon shows linear dependence as in- 
dicated by red lines. The shaded light blue regime represents 
the continuum of magnons in the bulk. 



width Ly is sufficiently large so that the coherent overlap 
between opposite edges can be ignored. 

Eliminating the field tp*^ with the help of Eq. (Is]) , the 
boundary conditions on the upper edge is simplified to 
the constraint on the field tpA solely, 



dy 



y/2-fJ Cb dt 



0. 



(7) 



Note that the above relation impose constraint on how 
the (imaginary) momentum in the transverse direction 
renormalizes the propagation of magnons on the upper 
edge. It can be shown that edge magnons on the upper 
boundary carry quantum number /S.Sz = — 1 with evanes- 
cent wave function (l>A{x,yjt) = 4>Ae{x,t)e°'y'^ , where 
ttj, > is the imaginary momentum along the transverse 
direction. Substituting the boundary constraint into the 
bulk K-G equation, the dimensionality is effectively re- 
duced to one. The resultant equation for edge magnon is 
the one-dimensional K-G equation. 



1^ 



dx'^ 



2 2 



(j)Ae{x,t) = 0. 



(8) 



The spin-wave velocity and the effective mass for the edge 
magnon are related to its bulk values. 



Cb 



rub 



1 



< 1. 



1 + 



(9) 



27= 



The above results are plotted in Fig. 2. Since the excita- 
tion gap nieC^ and the spin- wave velocity Ce are smaller, 
the dispersion for the edge magnon lies below the con- 
tinuum and remains sharp even when the interactions 
between magnons are included perturbatively. 

On the lower boundary, the edge magnons carry quan- 
tum number AS'z = 1 with wave function 4'*g(x,y,t) — 
4'*Bei-'^'^)^~°'^^- Following similar calculations, one can 
show that the edge magnon also satisfies the one- 
dimensional K-G equation with identical parameters. 



The similarity between the edge magnons on opposite 
edges calls for a Z2 symmetry argument. It turns out 
the discrete symmetry relating these evanescent modes 
is the parity symmetry Py in the transverse direction. 

Because the Neel state has staggered spin configura- 
tion, the operation of parity symmetry needs extra cau- 
tion. When reversing the y— axis, it is clear that the 
lattice coordinates transform as (y. A) — )■ (— y, B) and 
{y,B) -^ {—y,A). However, due to the staggered spin 
configuration in the Neel state, the parity transforma- 
tion also reverses the spin orientations and causes "charge 
conjugation" effectively. That is to say, the parity trans- 
formation turns particle-like excitations into the hole-like 
and vice versa. Therefore, the solution under Py trans- 
formation takes the form, 



i)A{x,y,t) 
b*Bix,y,t) 



4>*B{x,-y,t) 

<j>A{x,-y,t) 



(10) 



The above symmetry is exactly what happens in the 
boundary field theory for edge magnons. 

The boundary field theory supplemented with the sym- 
metry argument fully answers our puzzle. The ferromag- 
netic magnons satisfies the one-dimensional K-G equa- 
tion originated from its bulk counterpart with explicit 
relations. The edge magnons running on upper and lower 
boundaries carry opposite quantum numbers and are an- 
tiparticles to each other related by the Py parity symme- 
try. In fact, the whole field theory (including the bulk 
and the two edges) is fully relativistic and excitations al- 
ways appear in pairs as required. The confusion mainly 
arises from the asymmetry of the spatial wave functions 
for the edge magnons because their antiparticles locate 
on the opposite edges. In short, the single-branch ferro- 
magnetic edge magnon on one zigzag boundary is indeed 
an antiferromagnetic one with its antiparticle running on 
the distant opposite boundary. The linear dispersion of 
the edge magnon (with specific relation to the bulk dis- 
persion) now looks more than natural. 

General cases. One may wonder how robust the above 
conclusions are if the Neel order in the bulk is absent. To 
explore this possibility, we represent the spin operators 
by Schwinger bosons which preserve the SU{2) symme- 
try in the isotropy limit 7 = 1 explicitly. Since the cal- 
culations are pretty much the same as the HP bosons, 
we skip the technical details here. We arrive at the same 
conclusions except the parameters Q and A in the Harper 
equations are modified due to quantum fluctuations. The 
robustness of the relativistic boundary field theory should 
not be a big surprise because the field-theory description 
for Schwinger bosons is still relativistic. 

We also perform DMRG calculations to demonstrate 
that the ground state is very close to the Neel state. It 
is known that the ground state for Heisenberg model on 
finite bipartite lattice is a spin singlet. Figure 3 shows 
the magnetization profiles {Sz(j)) in the lowest-energy 



(a) zigzag nanoribbon 




(b) rotated-square nanoribbon 



O : positive S' 



»: negative S' 0:15" I =0-01 



to explore and reveal the true nature of these edge exci- 
tations in the future. 
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FIG. 3: Magnetization profiles (5*2 (r)) for (a) in zigzag 
nanoribbon with Ly — 2 and (b) in rotated-square nanorib- 
bons with Ly — 3/2. Light (green) and dark (red) circles 
represent positive and negative values of the spin polarization 
respectively, while the areas are proportional to the absolute 
values. Crosses represent the edge site ro for which the local 
Zeeman field — ^^^(ro) with h — O.OIJ is applied. 



state of J2r ^z{^) = for the Heisenberg model with lo- 
cal Zeeman field —hSz{j^o) applied to a center site rg of 
one edge. The result suggests that even the small local 
field h — O.OIJ turns the ground state to almost perfect 
Neel state. It is also interesting to notice that the edge 
magnon is not a privilege of honeycomb nanoribbon with 
zigzag edges. For rotated-square nanoribbon as shown in 
Fig. 3, we repeat the same calculations and find the pres- 
ence of edge magnon as described in Eq. |8j with different 
parameters, 
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Ce 

Cb 






<1, 



1 + ^ 



(11) 



where Cb = ^Qa and rub — \/7^ — \/{AQa?'). Therefore, 
the emergence of the edge magnon is related to the un- 
compensated lattice structure. 

Although the Neel order may not play an essential role 
to the robustness of the edge magnon, one shall be cau- 
tious to draw similar conclusions for graphene nanorib- 
bon with the same geometry. Since edge magnetism in 
graphene nanoribbon is itinerant in nature, it is not yet 
clear whether the edge magnon can still be described by 
the one-dimensional K-G equation derived here. How- 
ever, recent Monte Carlo simulations 19, demonstrate the 
presence of sharp spin-wave excitation with sharp spec- 
tral weight, in qualitative agreement with otir boundary 
field-theory description. It would be of vital importance 
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